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CHAPTER 1 


INTRODUCTION 


The Hidden Line Problem 

To establish the need for this work, the author feels it is important that 
the reader obtain a basic understanding of the hidden line problem; this 
introduction is devoted to providing such an understanding. 

Figure 1-1 depicts a simple box prior to hidden line removal. Clearly it is 
difficult to get an accurate understanding of the actual orientation of the box; 
is it oriented as shown in Figure 1-lb or as in Figure 1-lc? Hence, one can see 
the great value in a hidden line representation of a model. 

The hidden line problem has been handled by many different individuals 
(l|. each with some uniqueness in their approach. But, aside from the 
uniqueness of individual approaches there are several considerations that are 
common to all hidden line removal methods. A number of simple illustrations 
will aid in a discussion of these similarities, and thereby provide some general 
insights into the hidden line problem. 

In the wire frame model of the box (Figure 1-la), each of its sides 
represents a surface; each of the surfaces has the capability of hiding any lines 
which may lie behind them. So the problem becomes one of determining 
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which surfaces are io front of which lines. To make this determination, a 
number of common testing procedures are used. 

First, and perhaps most simple, is the min-ma.x test. In Figure 1-lc it is 
apparent that side A can in no way obstruct line 1. The min-m.ax test would 
determine this by comparing the minimum and maximum coordinate values of 
surface A to the minimum and maximum coordinate values of line 1. This is 
analogous to drawing boxes around the coordinate extremes of surface A and 
line 1 as shown in Figure l-2a. If the two boxes do not intersect there can be 
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no intersection of the surface and the line, and therefore, the surface cannot 
hide the line. If the opposite is iiue, then there is a possibility that the 
surface and the line intersect, and further testing must be done to make this 
determination. 

This further testing comes in the form of what is called a "surrounded ness 
test". Figure l-2b demonstrates the need for such a test since the min-max 
test alone cannot eliminate th? possibility that the surface and the line 
intersect. The surroundedness test takes point "X" on line 1 and determines 
whether it is surrounded by surface A or outside of surface A; if it is outside of 
surface A, the line and the surface do not intersect, and if the opposite is true, 
it has been conclusively determined that the surface and the line do intersect 
in the X-Y plane. 

When it has been determined that the two elements intersect, all that 
remains is to determine whether the surface is in front of the line or vice versa. 
This is accomplished by th» "depth test". The depth test compares the 
minimum and maximum depth values between surface A and line 1; if the 
Ei'nimum surface depth is greater than the maximum line depth, then the 
surface is behind the line and cannot hide it (see Figure l-2c). On the other 
hand, if the maximum surface depth is less than the minimum line depth, then 
the surface will hide all or a part of the line (see Figure l-2d). 

The preceding has been a very simplistic dbcussion of the basic concepts 
in bidden line removal. There are many sorting methods employed in the 




different hidden line algorithms which speed up the above discussed processes. 
Further, the situations where a surface only partly hides a line or where 
surfaces penetrate into each other have not been discussed, yet must be 
considered. For a more in depth discussion of the hidden line problem the 
reader is referred to ”A Characterization of 10 Hidden-Surface Algorithms” [1] 




The Need for a Fast Hidden Line Algorithm 
With Optional Contours 

Perhaps the best evidence of the need for a fast bidden line algorithm in 
the MOVIE. BW display package is found in the speed comparisons between 
the modihed JonesD algorithm and the MOME.BYU modihed Watkin’s 
algorithm discussed in Chapter 6. The fact is, the Watkin’s Algorithm is a 
scan conversion algorithm which was written to facilitate the generation of 
continuous tone images, and as such, carries with it all the overhead required 
to produce such pictures.(3] This overhead, coupled with the overhead required 
to allow enhanced continuous tone images (ie. shadows and transparencies) has 
made the MOME.B^IJ modified version of the Watkin's algorithm quite slow. 
Further, when contours are desired, the MOME.B^TJ modified Watkin’s 
algorithm goes through a very time consuming Gouraud interpolation scheme 
to generatt contour segments [3]. On the other hand, the vector nature of the 
JonesD algorithm allows the minimization of interpolation requirements in 
contour generation and this greatly speeds up the process. In light of the 
above, the need for a fast hidden line algorithm with accompanying contour 
option becomes clear. 


CHAPTER 2 


THE JONESD ALGORITHM 


At this point a discussion of the JonesD Algorithm is in order. This 
discussion will include reasons for the Algorithm’s selection, a description of 
how the algorithm works, and finally a discussion of the algorithm's 
limitations. To aid the reader in obtaining a general understanding of the 
algorithm, Figure 2-1 has been provided. The author suggests that the reader 
take a moment to study the flow chart, as doing so will greatly aid in 
understanding what follows; in fact, each time the reader encounters a 
flowchart through the course of this work, this same suggestion b 
recommended. 

JonesD was Selected 

There are four reasons that the JonesD algorithm was selected. The first 
and primary reason was the speed its author, Gary Jones, attributed to it. He 
did some quite elaborate testing of hb algorithm against other hidden line 
algorithms and showed that hb algorithm was faster in ever)’ case (4). Second, 
the JonesD algorithm was written specifically to handle the hidden line 
problem which aided in making its claimed speed believable. Third, because 
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the algorithm made provision for line elements, it was nicely suited for 
handling contour vectors. Finally, although the JonesD algorithm had some 
limitations, which will be discussed later, these limitations were considered 
either tolerable or surmountable. 
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How the JonesD Algorithm Works 

Preprocessing Requirements 

The JonesD algorithm requires signihcant preprocessing which results in 
more than worthwhile time savings downstream. The basic components of 
this preprocessing are (1) the hashing of redundant nodes, (2) the creation of a 
vector list from the model element connectivity or surface list and (3) the 
hashing of redundant edges. A more detailed explanation of these 
preprocessing methods and the reasons for their use follows. 

(1) Hashing Redundant Nodes 

The hashing algorithm employed in the elimination of redundant nodes 
takes the X, Y, and Z coordinate values of each node, and from these values 
generates a key into the hash table. The bash table internal numbers are then 
used as indices into an array that contains the node numbers for the 
coordinates. Now, each time that the same key is generated a redundant node 
will have been found. So rather than enter the redundant coordinate into the 
coordinate array, it is skipped, and the hashing internal number is then used 
to retrieve the node number previously assigned to the coordinate values. The 
node number retrieved is then used in the connectivity list for the element 
being processed. Once the redundant nodes have been eliminated the next 
preprocessing step, the creation of a vector list, is possible. 
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(2) Creating the Victor List 

The vector list referred to here is simply an array which contains pointers 
into the JonesD coordinate array for the endpoints of each of the vectors in 
the model; this list is required by the JonesD hidden line processor which will 
be discussed later in this chapter. The vector list is created by looping 
through the element connectivity array and creating a new two dimensional 
array containing the endpoints (pointers into PGRID) for each vector in the 
model. Figure 2-2 illustrates this process. While the vector list b being 
created it is advantageous to flag redundant vectors, since such flagging 
dramatically reduces the amount of computation that the JonesD algorithm 
must do. 

(3) Hashing Redundant Edges 

The same hashing algorithm used to eliminate redundant nodes is 
employed in the flagging of redundant edges. The vector endpoints are first 
arranged in ascending order after which a hashing key is generated from the 
two endpoint values. For each time the same key b generated, a redundant 
edge flag (JREDUN) b incremented by 1 and assigned to the vector being 
worked with. Consequently, the vectors with a redundant edge flag greater 
than 1 are ignored (Figure 2-2) in the JonesD hidden line processor. It should 
be pointed out that while the edges are being flagged for redundancy, they are 
also flagged as being edges of surface elements (AT^TE). This b done to 
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FIGURE 2-2 

LOADING THE VECTOR LIST 


distinguish surface vectors from strictly line element vectors such as contours; 
this again speeds the operation of the JonesD processor. 

For more information about the hashing algorithm used in this work the 
reader is referred to [5] and Appendix C. 




II 


JoncsD Hidden Line Preprocessing 

In addition to the preprocessing needed to get the data required by the 
JonesD algorithm, the actual processor itself does a significant amount of 
preprocessing. This preprocessing involves bucket sorting, vector preparation, 
and surface preparation. 

Bucket Sorting 

The JonesD Algorithm employs a number of processes that increase speed 
of operation. One of these processes is the sorting of both surface elements 
and surface vectors into buckets. These buckets are actually cells in an X-Y 
grid. The number of grids that the screen is divided into depends on the 
number of surface elements and surface vectors in the model. Further, the 
bounds for the grid are determined from the minimum and maximum X and 
Y transformed model coordinates; this results, hopefully, in a fairly even 
distribution of surface elements and surface vectors among the grid cells 
(buckets). Figure 2-3 visually depicts the preceding discussion. All of the 
bucket sorting information required by the JonesD processor is set up in the 




subroutine GVTHGD. 
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Figure 2-3 

Spatial Cells used in X-V Bucket Sorting 
From Reference [4] 


Vector Preparation 


Following is an explanation of the vector preparation that is performed in 
subroutine G\TH\P and diagramed in Figure 2-4. The process begins by 
loading the \TCTOR array with information about each of the vectors that 
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has not been flagged as redundant; this information includes endpoint 
coordinates, line equation information, and minimum and maximum 
coordinate values for each vector. Tue surface vectors are then mapped into a 
vector map (JX\T) according to the celb (buckets) that their endpoints begin 
and end in, along with any cells that the vector could potentially cross (Figure 
2-3). All of the vectors, surface or line elements, are them mapped into an 
additional array (J\T). Finally, the vectors in the surface vector map are 
sorted according to depth beginning with the vector closest to the eye of the 
observer and progressir.g back from there; a Shell sort is used to accomplish 
this. The reasons for this vector preparation will become evident in the 
description of the JonesD hidden line processor. 
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Surface Preparation 

A brief explanation of the surface preparation performed in subroutine 
G\THSP and diagramed in Figure 2-5 follows. This process begins by loading 
the XP, "^T, and ZP arrays with the coordinates of the 6rst element being 
processed. The maximum and minimum X, Y, and Z coordinate values arc 
then determined and loaded into the SITIF array. Once this is accomplished, 
the surface is mapped into the surface map (JX^’S) based on its minimum and 
maximum values in X and Y (Figure 2-3); these minimum and maximum 
values determine which cells (buckets) the surface might lie in. Next, the area 
of the surface is computed by taking cross products between nodes, and is 
then loaded into a location in the SlTtF array. During the area calculation, 
the component areas from each cross product are stored in the ,AR array for 
later use in the JonesD processor. The above steps are repeated for each 
element. The final step in surface preparation involves the sorting of the 
surfaces according to their depth; a Shell sort is used which sorts the surfaces 
from front to back according to their minimum Z coordinates. As with the 
vector preparation, the reasons for the surface preparation will become evident 
in the discussion of the JonesD hidden line processor. 
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For a more detailed discussion of surface and vector preparation, the 
reader is referred to the comments in subroutines G\THSP and G\Tff\T in 
Appendix A. 


The JonesD Hidden Line Processor 

The JonesD hidden line processor, which can be found in subroutine 
G\THID, can be divided into two basic parts. The first part computes vector 
intersection locations, and the second part computes the visibility of line 
segments and outputs them to the graphics display device. 

Vector Intersections 

The flow chart, Figure 2-6, visually illustrates the processes involved in 
the computation of .vector intersection locations, and the segmentation of 
vectors at such intersection points where required. As can be seen in Figure 
2-6, vectors are handled one by one starting at the beginning of the vector list 
(ATCLS) and proceeding through to the end of it. The main function of this 
section of the processor is to determine if any surface vectors lie in front of the 
vector currently being processed; if this is the case, the vector being processed 
must be divided at all points of X-Y intersection with surface vectors. The 
reason for this, as shown in Figure 2-7, is that every vector that crosses an 
element boundary must ne divided at that boundar}*; this insures that all 
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VECTOR INTERSECTION FLOWCHART 
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vector segments will be completely in front of, or completely behind surfaces. 










i: 


The vector jntorscctioD calculation bci;ins by locating all surface vectors 
that lie in the same bucket or buckets as the vector being processed (active 
vector), and which are in front of the active vector. This is easily 
accomplished because of the vector mapping and sorting, along with the Z 
coordinate minimum and maximum determination, that was done for each 
vector in subroutine G^T^^^T. A binary sea.’ch locates, from the sorted 
surface vector map (JX^T), the dividing point between surfaces that lie in 
front of, and behind, the minimum Z coordinate )f the active vector. This 
idea is shown graphically in Figure 2-8a. Once the appropriate surface vectors 
have been located, each one must fail two trivial rejection tests before an 
intersection location will be calculated. The first of these two tests is the 
min-max test that was discussed in the introduction to this work. If the 
bounding box around the surface vector being looked at, and the bounding box 
around the active veclor do not overlap, ther^ b no need for further testing of 
the current surface vector and the algorithm goes immediately to the next one. 
If, on the other hand, the min-max test reveals a possible intersection between 
the two vectors, a second rejection test b employed. Thb test compares the 
endpoints of the two vectors since the surface vector and the active vector 
could either be one and the same, or could be connected to each other. In 
either case there would be no point in computing the intersection between the 
two If the surface vector b not rejected by either of these two tests, 
calculations to determine a possible point of intersection between the two are 
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begun; these calculations are accomplished via the line equation data for each 
of the vectors that was created in subroutine G\TIf\T. As possible 
coordinate intersections are determined, they are checked against the 
minimum and maximum values belonging to the active vector. If these 
minimums and maximums lie outside the computed coordinate intersection 
location, the surface vector can be rejected (Figure 2-8b), Finally, if a point of 
intersection is located, the active vector b divided into two line segments at 
that point of intersection. Once the active vector has been divided into 
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segments, a Shell sort is once again employed to order the segment list from 
minimum to maximum, and the second logical division of the algorithm, 
segment visibility testing, begins. 

Segment Visibility Testing 

At the commencement of this discussion it is important to realize that the 
segment visibility testing process is handling the segments that were created 
by dividing up one vector in the vector intersection portion of the code. In 
other words, one vector at a time comes down the hidden computation 
pipeline, in which it is first divided into segments; each segment is then 
visibility tested and flagged as invisible if such is the case. If a segment is still 
visible it is immediately drawn on the graphics device, after which the hidden 
process is started over on the next vector in the vector list. 

.\t this point it will prove informative to follow one segment through the 
visibility testing process. First, the midpoint coordinates of the segment are 
computed from its endpoints; these midpoint values will be used for 
comparison purposes. This can safely be done because all vectors have been 
divided at locations where they intersect surface edges (Figure 2-7). This 
division insures that if the segment midpoint is behind a surface, the whole 
segment must be behind the surface and vice versa. Now, knowing the 
midpoint location, the surface bucket in which it resides is determined. A 
binary search is then used to find the point in the sorted surface list that 


divides the surfaies that fall in front of the segment from those that lie behind 
it (Figure 2-9). The segment is then checked against all of the surfaces that 
could lie in front of it, and therefore hide it. The following tests are used to 
accomplish this visibility check. First, as in the vector intersection portion of 
tie processor, the min-max test is used. Once again, if the minimum and 
maximum X and Y coordinate values of the surface and the vector segment do 
not overlap, the surface is trivially rejected since it could not possibly hide the 
segment. If the surface cannot be rejected solely on the basis of the min-max 
test, a more rehned containment test is done. Figure 2-10 illustrates the basics 
underlying this test; it is accomplished by computing the areas of all the 
triangles, created by joining each of the element vertices to the midpoint of 






21 







/ 

/ 



/ 

POINT A CONTAINED OY 
ELCNCNT APCA - APE A 
or SUODlVtOSO TRIANGLES 
* AREA or WHOLE ELEMENT 

POINT A OUTSIDE Or 
ELEMENT - AREA Or 
SUBOIVtOKO TRIANGLES 
» AREA or WHOLE ELEMENT 

FIGURE 2-iOA 

FIGURE 2 - 

iOB 


the lice segment, as shown in Figure 2-lOa. A ratio of this computed area 
with the area of the element, which was calf'ulated in the surface preparation 
routine GVTHSP, is then taken. If the point lies within the element this ratio 
will clearly equal 1.0 (assuming convex elements) as Figure 2-lOa shows. If 
however, the point lies outside the element, the ratio will be greater than 1.0 
as is clearly illustrated in Figure 2-lOb. Further it will have been determined 
\^hether or not this surface hides the line segment being processed. If it does 
not, the surface will be rejected, and the next potential hiding surface will be 
considered. If it is determined that the surface does surround the segment 
midpoint, an initial depth test is performed to see if the minimum Z 
coordinate of the surface is greater than the Z coordinate of the midpoint. If 
this is the case, the surface is rejected since the segment is in front of the 
surface and cannot be hidden by it. If the possibility still exists that the 





surface is io front of the line segment, the next test performed is one that 
determines if the segment is actually one of the surface edges. Figure 2-11 
illustrates how this is accomplished. If the midpoint is on the surface 
boundary, then the area of one of the triangles, used in the containment test 
and as shown in Figure 2-11, will be 0.0; this indicates that the segment is a 
part of the surface edge, and as such, cannot be hidden by the surface. If the 
surface makes it through all of the above tests, the plane equation for the 
surface is determined from the coordinate values of its first three nodes. The 
X and Y coordinates of the line segment midpoint are then substituted into 
the plane equation; this yields the Z depth of the plane at the X and Y 
coordinates of the line segment midpoint. The Z depth of the line segment is 
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then compared to the Z depth of the plane, and if the plane is found to be in 
front of the segment, the segment is flagged as invisible. Once all of the 
segments which result from a single vector have been visibility tested, the 
visible ones are drawn, and the endpoints of the vector are flagged for 
visibility. 

The above process is repeated for each vector in the vector list with the 
end result being an accurate hidden line representation of the model processed 
if four conditions are maintained. These conditions give rise to a discussion of 
the limitations of the JonesD hidden line algorithm which will be handled in 
the next chapter. 

Limitations of the JonesD Algorithm 

Since Gary Jones wrote his algorithm to process finite element models for 
the NASTFL\N plotting package, there were four considerations in hidden line 
processing that he did not need to consider. First, there was no need to be 
concerned with more than four sided elements since the majority of finite 
element problems deal with three and four sided elements. Second, 
penetrating polygons were not considered since finite element models require 
nodes at all intersection points. Third, the capability of handling concave 
elements was not considered because, again, such elements are not common in 
finite element models. Fourth and finally, consideration was only given to the 
handling of planar elements. 



CHAPTER 3 


MODIFICATION TO THE JONTSD ALGORITHM 


As briefly mentioned at the conclusion to Chapter 2, there were four 
limitations to the capabilities of the JonesD algorithm. These included (1) the 
inability to handle more than four sided elements, (2) the inabiiity to handle 
penetrating elements, (3) the inability to handle convex elements and (4) the 
restriction to only planar elements. One of these limitations has been 
overcome by modifying the JonesD algorithm; suggestions for ways to 
overcome the remaining three have also been considered. A discussion of both 
implemented modifications and suggested modifications for the future follows. 


N'Sided Element Capability 


In its original form, the JonesD algorithm handled four sided elements 
and line elements. Triangular elements were treated as degenerate four sided 
elements with the fourth node being set equal to the first node. As a result of 
this four sided limit, Gary Jones was able to hard code a lot of his program; 
this means, that in the surface preparation, for example, he could prepare all 
the surfaces as though they had four sides; he could compute areas, minimum 
and maximum c<x)rdinate values and plane equation information from closed 
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form equations. Figure 3*1 shows an example of the hard coded form of a 
section of the JonesD algorithm, along with its counterpart n-sided code. 

In order to allow n-sided capability, one basic change in the code was 
required. This change involved keeping track of how many sides there were in 
each element in the preprocessing section of the code. This knowledge allows 
the use of ”do loops” in area calculation, minimum and maximum coordinate 
determination and in visibility testing. For example, the area of an element b 
computed by taking cross products around the element nodes; this is easily 
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accomplished in a loop that goes from one, to the number of sides in the 
element. The same type of loops were employed in the determination of 
element maximum and minimum coordinate values and in line segment 
visibility testing. This generalization of the code required the creation of three 
new arrays to handle the X, Y, and Z coordinates of each element; 
additionally, one new array was required to keep track of the partial areas 
determined from the nodal cross products for each element. Finally, the 
connectivity array was expanded to handle n-sided elements as required by the 
user. 

It is the authors opinion that a detailed description of each step taken to 
make the n-sided enhancement would not provide any new insights to the 
reader. The fact is, the actual modiScation b not where the greatest difficulty 
arose; the greatest difficulty came in the authors efforts to understand how the 
JonesD algorithm worked. This difficulty came as a result of the sparse 
comments contained in the original code. 

Penetrating Eleznenta 

As mentioned above, the JonesD algorithm cannot accurately handle 
penetrating elements. Figure 3-2 demonstrates the problems that the 
algorithm has. Notice that at the intersection of the two spheres, the 
algorithm does not know where to terminate the penetrating lines The reason 
for this is that the points at which lines intersect surface planes are never 
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computed in the line intersection portion of the code. That is, intersections 
are only computed at the locations where lines intersect lines, and not where 
lines intersect surfaces. Consequently, if the midpoint of a line, that 
penetrates a surface, is in front of that surface, it will remain completely 
vbible. Conversely, if the midpoint of a line, that penetrates a surface, falls 
behind the surface, the entire line will be rendered invisible. 

The apparent solution to this problem would involve the computation of 
the points at which penetrating lines intersect planes. This could be 
accomplished using the line equation of the "'penetrating line" and the plane 
equation of the element it penetrates. This would slow the algorithm down 
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considerahl) , since a line in a bucket would have* to be tested against ever)’ 
surface whose minimum and maximum coordinate values were in common 
with the minimum and maximum coordinate values of the line. Further, this 
would only result in the termination of the penetrating lines where they 
intersected surfaces. The lines of intersection still would not be drawn. 

At any rate the implementation of this capability is left to another 
investigation. 


Concave Elements 

As Figure 3*3 illustrates, the JonesD algorithm does not handle concave 
elements with any degree of reliability. This limitation results from the area 
ratio method, used in determining if a line segment's midpoint is contained 
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within a surface element. In Figure 3-3 the sum of Ihe areas of the triangles, 
created by drawing lines from point A to each of the element vertices, is 
clearly greater than the area of the element. Consequently, the containment 
test would conclude that the point "A” lies outside of the element and, 
therefore, must be visible. This is of course not the case and results in the 
hidden line remaining visible. 

It is the author's opinion that this problem could be handled with 
minimal time expense since convex elements are not used very often in 
modeling. The solution would involve the taking of dot products between 
element edges which could be accomplished at the time that element areas are 
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FIGURE 3-4 
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calculated in subroutine G\THSP. As an element was traversed, the dot 
products of the successive edges could then be compared, and if the sign of the 
product changed from one set of edges to the next, one would know that a 
concave element had been encoui.tered. With this knowledge, the element 
could then be divided into two hopefully convex elements as shown in Figure 
3-4, and the surface preparation could be done on the two new and now 
convex elements. If division of the element resulted in a concave and a convex 
element, then the remaining concave element could be divided in two. This 
process would continue until only convex elements remained. From this point 
the hidden line process would proceed as usual and an accurate hidden line 
view obtained. 

Again, as much as the author would like to undertake this modiGcation, 
time requires that it be left to the future. 

Warped Elements 

Clearly, the JonesD algorithm would have limited success in processing 
warped elements since it relies on the plane equation as the ultimate test for 
line segment vbibility. The solution to this problem would be very time 
expensive since it would require that all elements be tested to find those that 
are warped. Such testing could be accomplished by comparing the normals 
around an element. If the normals varied, a warped element would have been 
found, and the element would have to be divided into triangles. Each triangle 
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could then be accurately processed by the modified JonesD algorithm. It is 
interesting to note, hovsever, that as long as the elements are not too severely 
warped, JonesD does a very acceptable job (Figure 3-5). 

The implementation of a contour algorithm, as mentioned above, will be 
the subject of the next chapter. 



Figure 3-5 
Mt. St. Helens 


CHAPTER 4 


CONTOUR IMPLEMENTATION 


The implementation of the contour routine was one of the most 
interesting and rewarding parts of this work. The algorithm was conceived by 
an unnamed employee of the John Deer Company with information that was 
obtained from L. J. Segerlind’s book ”Applied Finite Element Analysis" |6|. It 
is very simple, yet produces high quality results. However, to obtain these 
high quality results, a number of modiGcations and enhancements had to be 
made to the basic algorithm. These modifications can be adequately discussed 
in five divisions which are, (1) the justification for computing the contours 
prior to entering the hidden routine, (2) the contour problem, (3) the contour 
generation algorithm, (4) contour sorting and (5) contour hidden processing, 
output, and labeling. 

Justification for Contour Preprocessing 

By glancing at Figure 4-1, the reader will quickly be able to see where the 
contour package fits into the overall hidden line processing system. Now, 
realizing that many of the contours that are created on a 3-D model are going 
to be hidden, one mi^ ,t think that generating contour segments prior to 


h 


33 


31 



riouRe 4~! 

CONTOURS IN HIDDEN FLOWCHART 


hidden line removal is wasteful. However, there are three good reasons for 
generating contours before hidden processing. First, and most convincing, is 
that there b no accurate way to create contours given the strictly vector 
output from the hidden routine. The problem b that once returned from 
hidden, all there is to work with are coordinates in 2-D space. There b no 
way to tell, as shown in Figure 4-2, which surface lies in front of which, and 
therefore, there is no way of knowing how to connect equal function locations 
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to obtain accurate contour lines. Second, the riumber of contour segments 
generated on a typical model is usually quite small relative to the number of 
lines that combine to make up that model. This normally small number of 
contour segments is handled quickly by the JonesD hidden algorithm. This is 
because the algorithm knows that these segments can only be hidden and that 
they cannot hide anything. Finally, prior to hidden processing, the element 
connectivity is intact; this makes the accurate generation and proper 
connection of the contour segments possible. 

The Contour Problem 

Before going into a detailed discussion of the contour algorithm, a basic 
description of the contour problem is in order. Contours are created by 
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interpcdaling between the function values assigned to elemen* coordinate?, as 
diagramed in Figure 4-3. T^ese function values may represent any type of 
function, and are assigned to the the element nodes according to some type of 
preprocessing operation. The preprocessor might be a finite element package 
for example. If one wishes to ind the contour line at a function value of 15, in 
Figure 4-3, interpolation n» ?ds to be performed between adjacent element 
nodes to find all of the points along the clement edges that have an 
interpolated value of 15; these locations are designated with an ”X’. Once 
this is accompl’shc , ‘♦>e two ”X” points are connected producing a contour 
segment. 



FIGURE 4-3 
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The Contour Generation Algorithm 

The flowchart in Figure 4-4 diagrams the basic processes in the contour 
generation algorithm as implemented in conjunction with the JonesD hidden 
line algorithm. The initial portion of the algorithm divides each of the 
elements into triangles, by creating a center node which is the average of esch 
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of the element nodes. The average function value at the center node is also 
computed by averaging the function values at each of the element nodes. This 
process is done to transform the element data into the triangular form required 
by the contour routine. The main reason for using triangular elements, as a 
basis for computing contours, is that triangular elements possess a unique 
interpolation property. This is demonstrated in Figure 4-3 where it is seen 
that if a contour value is present along one edge of a triangle, it will also be 
present on one and only one of the other edges of that triangle. Conversely, if 
a contour value is not present on two sides of a triangle, it cannot exist on the 
third side. Because of this property one can be assured that contour lines will 
never cross within triangular elements. Further, one can be assured that 
contour lines will never cross over a series of elements if each of those elements 
has been subdivided into a triangles for contour generation purposes. This 
idea can be visually understood by viewing Figijre 4-5. 

The contour algorithm next begins to process each of the triangles that 
' make up the first element. One by one each triangle is sent through the 

interpolation portion of the code. In this portion, contour segments, 
representing all of the contour levels that have been user requested and that 
lie within the nodal function values for that triangle, are generated. The 
process by which this is accomplished is as follows. 

The function value at the first contour level is extracted from memory 
and compared to the function values at each of the nodes thut make up the 
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triaDguIar element being operated on. If the 6rst contour value lies outside 
the bounds of nodal function values, the element is in.n.ediately tested against 
the next contour value. If, however, it is found that the contour value docs lie 
along the first edge of the triangle, then a contour segment endpoint, on the 
edge, is interpolated from the edge endpoint function values and X, Y, and Z 
coordinates. The next edge of the triangle is then tested for contour value 
intersection, and if it is found, the other endpoint of the contour segment is 
interpolated along the second edge. The contour segment endpoints are then 
immediately stored into memory. If it is discovered that the contour value 
does not lie along the second edge, then a contour segment endpoint is 
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immediately interpolated along the third edge, and the contour segment 
endpoints are stored into memory. This is done since it is known that such 
an endpoint must exist because of the interpolating nature of the triangle. 
This process is repeated for each contour value until all of the contour 
segments on the triangle have been computed 

Succeeding triangles in the element are then processed in a similar 
manner until all of the triangles making up that element are completed. The 
same process is then repeated for the next element, and so on until all of the 
elements have been processed. It should be pointed out, that while the 
contour segments are being created, the number at each different level, as well 
as the locations of their endpoints, b tracked. Because of this it b possible to 
sort the presently random contour segments into nicely ordered contour loops. 

Contour Ordering ‘ 

Contour ordering b required for two reasons - first for accurate labeling 
and second so that curves can be accurately fit through contour loops if 
desired. When the contour ordering routine G\THCO b called, it b passed 
contour segment coordinate information sorted by contour level; the 
information is not, however, ordered within those leveb. It b also passed the 
number of contour segments in each level. With thb information a systematic 
process of sorting the contour segments into ordered loops b started. This 
ordering procedure takes advantage of the fact that adjacent contour segments 
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uill have like coordinate endpoints. Consequently, a search starting with the 
first contour segment at a particular contour level is begun in an attempt to 
6nd another one, at the same level, with a like endpoint. This ordering 
procedure diagramed in Figure 4-4, flows as follows. 

(1) The endpoints of the first contour segment, in the segment storage 
array, are stored in holding locations. The first contour segment will be 
referred to as the ”initial segment" hereafter. 

(2) The #1 endpoint (see Figure 4-6) is initially tested against the 
endpoints of both ends of all the segments in its same contour level. If it 
is found that one of the other segments, at this level, has an endpoint 
that matches that of the #1 end of the "initial segment", a logical is set 
to true. This logical indicates that the "initial segment" has connecting 
vectors at both of its ends. 
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(3) The #2 endpoint of the ”initial segment" is then compared to both 
ends of the other segments in its contour level. When a vector with a 
matching endpoint is found, the endpoint storage locations of this vector 
and those of the second vector in the list are swapped. The sorting then 
proceeds starting with the #2 end of the vector that newly occupies the 
number two location in the storage area; the #2 end of this vector is now 
compared to the remaining vectors in the list at this contour level until a 
connecting vector is found. Figure 4-7 will clarify this sorting process. 

(4) When the end of the contour loop, which progressed from the #2 end 
of the "initial vector", is encountered, the array order of the contour 
segments is reversed if the logical indicating that there were contour 
vectors proceeding from both ends of the "initial vector" is set to true. 

The process describ'd in part (3) then continues until the last contour 
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vector at this level has been sorted or until no more connecting vectors 
can be found. If, however, there is no vector that connects to the #1 
end of the ’’initial vector", and there are still more segments at this level 
when the end of the loop being processed has been reached, then two 
possibilities exist. The 6rst one b that the "initial vector" happened to 
be at one end of the contour loop which does not pose any problem. The 
second b that the contour loop closed on itself in which case there is no 
point in reversing the order of the sorted segments; this excess sorting is 
avoided by storing the #1 end coordinates of the "initial segment", and 
then checking to see if they are the same as the #2 end coordinates of 
the last segment found in the loop being processed. If they are the same 
then order reversal b bypassed. 

(5) If it happens that the end of the 6rst contour loop at this level is 
encountered before all of the vectors have been dealt with, then there is 
clearly a separate contour loop at this level. Thb is easily dealt with by 
making the next vector, in the yet unsorted list, the new "initial vector". 
When this is done the process continues as described above until all the 
contour segments at a contour level have been processed. After this, the 
successive contour levels are processed until all of the contour segments 
have been ordered. 


During th** ordering process counters keep track of how many separate 
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coDtour loops there arc at oarh contour level, and hov> many contour se;;rnents 
are contained within each loop. This information is used in interfacing the 
contour sub-system to the hidden line sub-system. 

Contour-Hidden Line Interface 

There are two parts to the contour-bidden interface; one involves loading 
the contour information into the hidden algorithm, and the other involves 
correctly labeling the visible contour segments after hidden line processing. It 
turns out that once the contour segments have been ordered, the problem of 
loading them into the hidden line subsystem is quite simple. All that is 
required is that the contour vectors be added to the model vector array and 
that their coordinates be added to the model coordinate array. However, 
rather than load duplicate coordinates fr»^m the like ends of connecting 
vectors, the knowledge of how many segments are in each loop and of how 
many loops are at each level is used to allow the entering of coordinates only 
once. This loading operation is performed in subroutine GNPHCL. Now that 
the contour vector information is in the model vector and coordinate arrays, 
the hidden algorithm processes them as though they are a part of the model. 

The second part of the contour-hidden interface deals with keeping track 
of where contour loops become invisible due to hidden line processing and 
where they once again become visible. This operation is done at the time that 
line visibility is determined in the hidden algorithm. Once again, counters are 
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employed to keep track of how many separate contour loops there are after 
hidden line processing. The number of segments in each of these loops is also 
counted. First, this allows curves to be 6t through contour loops. Second, it 
rllows for efficient labeling of contours, since each time a contour loop becomes 
visible, it is labeled. That is, each time a contour loop becomes visible, the 
screen coordinates of the 6rst segment, along with the contour level number 
for the segment, are sent to the MOME.B^T LABELS subroutine. As a result, 
ever)’ visible loop in ever)’ contour level is labeled. Figure 4-8 shows the 
results of this proces.ses. 



Figure 4-8 
Contoured Cube 


All that has been discussed from the beginning of this work to this point 
constitutes the main body of new innovations in this thesis. The remaining 
two chapters deal with the interfacing of this work to the MOME.B'^T* 
Display package and with epu testing of the algorithm to check its actual 


CHAPTER B 


INTERFACE TO MOVIE.BYU DISPLAY 


The challenge in interfacing the modified JonesD algorithm to the 
MOME.B^XI Display package was in pulling out the polygonal data at the 
correct location and getting it into the correct form. This required the 
addition of eight new routines that closely parallel eight that are already 
present in MOME.BW Display. It also required the addition of one 
completely new routine to order the edges of clipped polygons. An explanation 
of the interface and the reasons for the methods used to accomplish the task 
follows. 

Data Required by the JonesD Algorithm 

The first thing that was considered in creating the interface was the data 
that the modified JonesD algorithm required. Since the JonesD algorithm’s 
only function is to determine hidden line removal, it was clear that all 
polygonal data had to be picked up after transformation and clipping 
operations had been performed. Therefore, the three elements of data required 
by the JonesD algorithm which are, (1) the element coordinate array, (2) the 
element connectivity array along with the number of sides in each element and 
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(3) the function or contour information, ucro all created or obtained after 
MOME.B^TI Display routines had performed all the necessary tr. iSformation 
and clipping operations. 

Modifleatior > to MOVIE.BYU Display 

The first interface to MOVIE. Display was the ’quick fix’ type where 
one obtains the information without regard for what has taken place in the 
code prior to the time it is picked up. This type of an interface, although it 
may work, usually accomplishes three things. First, it makes the existing code 
less understandable; second, it makes the program run slower since invariably 
a lot computation, not needed by the new- subsystem, has been done; and 
finally, it removes the desirable feature of modularization. As a result of these 
considerations, it was determined that a more complete interface would be 
attempted. This resulted in some close scrutiny of all the routines that 
preceded the location at which the required data was obtained. Because of 
this scrutiny, eight new routines were created which represented only small 
part!; of eight MO\TE.BVJ routines. These new routines, with thoir 
MOME.B'^X’ counterparts, are diagramed in Figure 5-1 and are contained in 
Appendix D. The reduction in routine size came mainly from the deletion of 
all of the color, shadow, and transparency information contained in the 
MOME.B'^X' routines, and resulted in a clean, efficient, and modular interface. 
This interface allows the transformations, clipping, poorman’s hidden line 
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removal, shrink option, and contour labeling currently available in the 
MOME.B'^X' Display package. It does not presently allow node numbering, 
element numbering, and the dotted line capability, although it could be 
modified to handle such. 
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Obtaining the Clipped and Transformed Date 

The clipped and transformed data is obtained in the new C\TIIPS 
routine which is a counterpart to the MOME.B^Xi POLSNT routine. This 
routine calls the G\PHCL routine ( counterpart to CLIP) then loads the 
clipped coordinates, connectivity, and function information into holding 
arrays. This operation proceeds on an element by element basis with all of the 
edges of one element being successively clipped. The edge information is then 
stored before the next element is processed. Because the polygon edge clipper 
in MOME.B'^X^ does not return the polygon edges in consecutive order, they 
are ordered (G\THLD) before they are loaded into the holding arrays. Once 
all of the element edges in the model have been transformed and clipped and 
the holding arrays have been correspondingly Slled, the hashing preprocessor 
to the modified JonesD algorithm is called and operates as was discussed in 
Chapter 2. From the" preprocessor, the model information goes to the modified 
Jone.^D bidden line processor which calls the present MO\TE.B'\X' PLTLIN 
routine to output the visible model segments. If contours are enabled, the 
MOME.B^X L-^BELS routine is called to label the contours. 

Now that the modified JonesD algorithm with all of its appendages 
(Figure 5-1) has been discussed, all that remains b to test the algorithm to see 
if it is really as fast as its author has claimed \i to be. 


CHAPTER 0 


SPEED COMPARISONS AND CONCLUSIONS 


.AJt hough the speeds recorded by the JonesD algorithm were not as good 
as those indicated by its author (4], it did perform ver>’ well in a number of 
cases. It was compared only to the MOME.B^Tj Watkin's algorithm. This 
comparison included the time necessary to transform and clip all models, and 
the time to perform the actual hidden computations and output the visible 
segments to the graphics device. If contours were computed, it also included 
the time necessary to generate them. All models were tested in exactly the 
same manner for both the Watkin’s and JonesD algorithms; that is, regardless 
of the algorithm, the same command 61es were used in generating hidden line 
pictures. The command 61e generated pictures are included in Figures 6-2 to 
6-11 at the end of this chapter. 

The items of interest in making speed comparisons were, (1) how would 
Watkin’s and JonesD compare as the number of model elements was increased, 
(2) how would the two algorithms compare in the area of contour generation, 
and (3) how would they compare if a model had high concentrations of 
elements in localized areas. 


51 


Inrrca.'if Number of Elomcnts 


To determine the effect of successively increasing the number of elements 
on algorithm speed, a prismatic hexahedron was selected. The number of 
elements per hexahedron face was systematically increased from 4 to 2S9 (see 
Figure &-2 and Figure 6-3). The results of this testing are contained in Figure 
6-1 and indicate two things. First, increasing the number of elements between 
the ranges of 24 and 864 only slightly affects the speed of the JonesD 
algorithm. The vectors processed per cpu second bares this out. However, 
when the number of elements exceeded 1000, the speed (see Figure 6-1) of the 
JonesD algorithm quickly dropped off. On the other hand, the Watkin's 
algorithm started out verj- slowly, and then picked up speed as the number of 
elements was increased. Consequently, the speed advantage of the JonesD 
algorithm was lost in the 900 to 1000 element range. 

Contour Generation 

As anticipated, the modified JonesD package is much faster at generating 
and displaying contours than the Gouraud interpolation method used in the 
Watkin’s algorithm. For small models, the contour generation and hidden 
representation capability is three to four times faster in the modified JonesD 
(Figure 6-1). Figures 6-4 and 4-8 are examples of these small models. As 
models increase in size, the speed of the JonesD package slows down (Figure 
6-5). In the area of generating contours on a flat grid as might be the case in 
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a two dimensional 6nite element analysis, the JonosD algorithm is four to five 
times faster. This time savings is demonstrated in the Nfounl St. Helens 
model (Figure 5-6) where contours are generated on the flat grid which, when 
warped, becomes Mount St. Helens (Figure 3-5). Two actual finite element 
models, Figures 6-7 and 6-8 demonstrate similar time savings. In addition to 
the time savings, a comparison of Figures 6-8 and 6-9 shows no great difference 
between the JonesD and MOME.B'^’U contours. In fact, the only difference is 
that the JonesD contours, because of there vector nature, are rougher than the 
MOME.B^X^ contours; this problem diminishes as element size decreases. 

Element Concentration 

The Arne’s Plane model (Figure 6-10) was used to test the effect of 
concentrating large numbers of elements in localized areas of the viewing 
window. The wheels of the plane are each made of a large number of small 
elements which satisfy the element concentration requirement. As Figure 6-1 
shows, the Watkin’s algorithm wa.s slightly faster than the JonesD algorithm, 
even though there were only 587 elements in the model. This slowing of the 
JonesD algorithm is a result of a large number of surfaces, and vectors, being 
mapped into the same X-Y cell for bucket sorting purposes (Figure 2-3). 
When this happens, the sorting processes used to make the JonesD algorithm 
faster become inefllcient since the sort depth in a bucket b very deep. Figure 
6-11 illustrates this same problem. 


Conclusions 


As a result of this work two conclusions become evident. First, the 
modified JoncsD algorithm provides rapid hidden tine processing for relatively 
small models (ie. 1000 elements or less). This speed is however dependent on 
the localized density of elements. If elements are spaced evenly throughout 
the viewing window, the algorithm is fast, and if there are concentrations of 
elements, the algorithm slows down. Second, the algorithm, in conjunction 
with the contour generation algorithm, provides very rapid contouring 
capability. Both of these conclusions present the modified JonesD package as 
ideal for moderate sized models and in particular, finite element models. 
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A FAST IIIDDEN LINE ALGORITHM 


WITH CONTOLTl OPTION 


Ronald B. Thuc 

Dppartmenl of Civil Enginffring 
M.S. Dfgrff, Drcembfr 1984 


ABSTRACT 


Tbfre is an on going desirf in the computer graphics community to develop 
faster hidden line algorithms that can give accurate hidden line representations of 
geometric models. This thesis presents the JonesD algorithm, modihed to allow the 
processing of N-sided elements and implemented in conjunction with a 3>D contour 
generation algorithm. The total hidden line and contour subsystem is implemented in 
the MOV’IE.BVU Display package, and is compared to the subsystems already 
existing in the MO\'IE.BYU package. The comparison reveals that the modified 
JonesD hidden line and contour subsystem yields substantial processing time savings, 
when processing moderate sized models comprised of 1000 elements or less. There 
are, however, some limitations to the modified JonesD subsystem. 
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